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We  review  the  recent  theoretical  study  of  the  effect  of  nonequilibri¬ 
um  phonons  on  hot-carrier  relaxation  and  transport.  In  a  quantum  well,  the 
proper  treatment  of  the  electron-phonon  coupling  between  clectroas  confined 
to  two  dimensions  (2-D)  by  phonons  traveling  freely  in  three  dimensions  (.J- 
D)  requires  special  care  because  phonon  heating  produces  a  bottleneck  in  the 
rate  of  transfer  of  energy  from  the  carriers  to  the  phonons.  Because  the  car¬ 
riers  interact  with  phonons  primarily  when  the  latter  are  close  to  the  quantum 
well,  the  latter  should  be  described,  not  by  plane  waves,  but  by  packets 
adapted  to  the  shape  of  the  carrier  confinement.  A  quasi-equilibrium  tech¬ 
nique  that  retains  off-diagonal  elements  in  the  phonon  wave-vector  permits 
an  unrestricted  treatment  of  the  density  operator  equation.  That  in  turn  leads 
to  a  choice  of  wave  packet  that  comes  from  solving  the  integrodifferential 
equations  rather  than  by  imposition.  Moreover,  if  the  carrier  distribution  is 
assumed  in  quasi-equilibrium  with  a  given  drift  and  temperature,  the  coupled 
partial  differential  equations  are  reduced  to  coupled  ordinary  differential 
equations  that  can  be  solved  with  modest  computer  power.  Comparison  with 
experimental  results  for  steady  flow  of  energy  from  carriers  to  phonons,  and 
for  time-dependent  relaxation  yields  quantitative  agreement. 
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1.  THE  PROBLEM 

The  pnxiuciion  ot  tnicrcielectronic  devices  by  molecular  beam  epita-xy.  and  the  high 
mobility  ol  ctuners  ptirticularly  in  modulated  heterostructures  created  the  importance  ot  this 
area  ot'  study.  Moreover,  tor  transport  in  small  high  mobility  systems,  inixierate  voltages  can 
lead  to  strong  tields  and  non-linear  effects.  An  excellent  review  of  two-dimensional  transpon. 
with  an  extensive  list  of  references  has  been  given  by  Ando.  Fowler,  and  .Stern' 

There  already  was  expenmenial  evidence  by  Shah  ei  and  Ryan  ci  at'  that  the 
rate  of  energy  transfer  from  electrons  to  phonons  was  an  order  of  magnitude  less  than  pertur¬ 
bation  theory  would  yield. ^  ^  On  the  other  hand,  controversial  experimental  results  were  also 
reported  by  A  ang  vt  al.^  Several  causes  that  may  contribute  to  this  reduction  of  the  energy- 
loss  rate  were  considered.  The  tirst  one  is  ihe  effect  of  reduced  dimensionality.  Theory' 
showed  merely  insignificant  dependence  of  the  energy-loss  rate  on  dimensionality,  when  the 
predominant  loss  is  by  emission  of  longitudinal-optical  (LO)  phonons.  The  second  is  screen¬ 
ing  of  electron-phonon  interaction  by  free  earners.  Screening  has  been  appraised  by  Riihle 
et  at  to  play  only  a  minor  role  up  to  electron  densities  of  n  =  lO’^  cm  \ 

It  was  clearly  perceived  by  the  expenmemalists.  and  a  qualitative  theory  was 
developed  by  Price,'”  that  the  inability  of  the  longitudinal  optical  (LO)  phonon  system  to  dis¬ 
sipate  heat  fast  enough  was  creating  a  bottleneck.  A  reabsorption  of  phonon  energy  by  elec¬ 
trons  is  presumed  to  decelerate  carrier  cooling.  In  fact,  this  bottleneck  effect  has  been  dis¬ 
cussed  by  Potz  and  Kocevar"  in  the  3D  case.  In  the  2-D  case.  Price's  theory  ts  qualitative, 
because  he  was  forced  to  introduce  an  ad  hoc  parameter  -  the  number  of  phonons  that  interact 
with  an  electron;  Price  recognized  the  need,  and  called  for  a  more  rigorous  treatment. 

2.  OUR  RESOLUTION 

During  that  ttme,  I  (M.  Lax)  perceived  a  need  for  understanding  electron  (and  hole) 
transport  in  quantum  wells  and  heterostructures  and  hired  two  research  assixiates,  W.  Cai  and 
M.  C.  Marchetti,  to  work  in  this  area.  Cai  was  already  an  expert  in  semiconductor  physics 
and  Marchetti  an  expert  on  transport  in  liquids.  My  role  was  advice  and  cnticism. 

We  concluded'*  '^  that  the  problem  was  a  general  one:  how  should  electrons,  whose 
!r;insport  is  confined  to  two  dimensions,  interact  with  phonons  that  can  propagate  freely  m 
three  dimensions.  .More  specifically,  since  the  electrons  only  interact  with  phonons  when  they 
.ire  in  the  vicinity  ot  the  ijuantum  well,  a  plane  wave  descnption  tor  the  phonons  ic  incon¬ 
venient.  One  possibility  is  to  retain  the  plane  wave  descnption.  but  quantize  the  phonons  iin 
the  :  direction)  over  a  thickness  /,  comparable  to  the  well  width.  But  such  a  treatment  would 
he  equivalent  to  Pnee's  with  /,  as  the  arbitrary  parameter. 
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Cai  proposed  a  resolution  of  this  problem  by  using  a  basis  set  for  the  z  direction  (nor¬ 
mal  to  the  well  walls)  that  consists  in  a  Gaussian  times  a  set  of  Hermite  polynomials.  But  it 
was  not  known  how  manv  terms  were  needed.  The  use  of  a  single  term,  the  Gaussian,  has  as 
an  arbitrary  parameter,  the  width  of  the  Gaussian. 

2.1.  Quasi-Equilibrium 

After  reviewing  the  work.  1  suggested  that  the  shape  of  the  phonon  ‘packet  ’  should 
come  out  of  the  problem,  not  be  imposed.  Marchetti  then  suggested  the  use  ot  a  quasi- 
equilibrium  procedure  of  the  sort  introduced  by  Bogolyubov  in  dealing  with  classical  liquid 
transport  and  by  Zwanzig'*  '"’  in  a  vanety  of  problems.  .See  also  Zubarev'^  and  Peletminskii 
.ind  Yatsenko''  The  basic  idea  is  the  quasi-equilibrium  assumption  that  the  vanables  ot  a 
problem  can  be  divided  into  slow  variables  and  fast  ones.  The  last  variables  are  assumed  to 
be  in  equilibrium  with  the  current  values  of  the  slow  variables.  When  ihe  fast  variables  are 
inserted  into  the  equations  for  the  slow  ones,  we  get  an  effective  set  ot  equations  tor  the  slow 
variables. 

2.2.  The  Choice  of  Slow  Variables 

The  success  of  such  a  procedure  clearly  depends  on  the  appropriate  choice  of  slow 
vanables.  The  electronic  variables  describe  the  occupancy  of  a  state  of  transverse 

momentum  k  in  the  nth  subband.  The  average  of  this  set  of  variables 

/«kU)  =  <ank‘'nk> 

is  the  familiar  distribution  function  for  these  carriers.  These  variables  must  clearly  be 
included  in  the  .slow  set. 

The  phonons  are  described  by  the  three-dimensional  wave  vector 

Q  =  (q.</z) 

associated  with  a  plane  wave  representation  exp(tQ  R).  The  variables  to  be  used  in  this  case 
are  with  average  value. 

Marchetti  made  the  crucial  proposal  to  retain  the  off-diagonal  elements  in  Eq.  (2..^). 
.Although  the  work  stans  in  the  plane-wave  representation,  by  allowing  off-diagonal  elements 
with  respect  to  </,.  we  have  prepared  the  way  for  an  eventual  transformation  to  packets  whose 
shape  is  as  yet  unknown. 

2.J.  Description  of  Ihe  Hamiltonian 

The  Hamiltonian  H  consists  in  an  electron  part  a  phonon  pan  an  electron- 
phonon  interaction  \  ,.p  and  a  phonon-phonon  interaction  \  ,,p  w  ith 
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He  =  'LEnk  tink  (2.4) 

where  is  the  energy  associated  with  the  state 

7)  = /)'‘'^Cn(-)exp(ikT)  ,  (2.5) 

associated  with  the  transverse  wave-vector  k,  and  quantum  well  state  n.  Here  C„(r)  is  the 
n-th  quantum  state  in  the  well.  The  unperturbed  phonon  part  is 

Hp  =  hQ^Q  (2.6) 

Q 

Q  =  (q.  iiz) 

The  electron-phonon  interaction  is  given  by  the  Frdhlich  interaction 

i’ep  =  jdRldR'e  n{R)  V'  P(R')|  ,  (2.7) 

namely,  the  Coulomb  interaction  between  the  electron  charge  density  en  and  the  charge  den¬ 
sity  -  V  P,  where  P  is  the  phonon  induced  polarization,  so  that 

^ep  ~  I  S  S  S  1  ^k',  k+<|  (j/i'nlq*  ~  ^k'.  k-q  CJ';,(q»  <7z)  f  ^n'k'^nk  > 

n.  k  n'.  k'  Q  I  J 


(2.8) 


where  the  matrix  element  G  takes  the  form 


G„  „{q,  qz)  =  j  dz  (,*'{z)  exp(  iq^z)  ^(z).  (2.9) 

Q'^L  - 

In  Eq.  (2.8),  a  is  the  Frdhlich  interaction  constant'* 

a  =  |2jte^Hcot(l/e„  -  l/En  (2.10) 


2.4.  Form  of  the  Equations  of  Motion 

One  can  derive  the  kinetic  equations  for  a  set  of  macroscopic  observables,  {Yi(t)}  from 
the  quantum-mechanical  Liouville  equation.  Here  we  choose  the  following  observables: 
l7i(0}  =  {  fnk(‘').  ^ •}  (see  Eqs.  (2.1)  and  (2.3)).  If  the  Liouville  equation  is  solved 

to  the  lowest  order,  the  rate  of  change  of  any  function  of  7(r).  denoted  by  F{y{t)),  is  given  by 


dFiyit)) 

rU 


^Trl  Hn(t).  F{^t))  lpo(t)-H 
h 


i 

Ji 


lim  j  dxe^TrW  V(T),  I  V'.Ffy)) )  po(r) }  , 


where  V'(t)  is  the  coupling  perturbation  in  the  interaction  representation: 


(2.11) 
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\'(T)  =  expl  Ui  fi)(He+Hp)x\  V  expl  (-i  riMH^+Hpn]  V 


12.12) 


and  po(f)  the  unperturbed  density  matrix.  Here,  neither  the  electron  nor  the  phonon  system  is 
assumed  to  be  even  close  to  equilibrium.  F-rom  Fiq.  (2.1 1).  we  obtain  coupled  integrodilYeren- 
lial  equations  of  the  form 


^fnk 

IT 


SLX  S  ■'^\fnk  -fn  k-  .  /  )1 

<1  <f!>h'nk' 


(2.1.^) 


=  S  Z  l.fnk  I  >-/n  k- 1  I  r'l'iir  terms  | 

i,,p  nkn'k’.;/’ 


where  a  typical  term  is  given  by 


5k',  k-q 


Gn'ntq.  t/,)  G;i'„(q.  q,"  ) 

e  -  l|£„k  -  £n-k-  -  fiUiL] 


q/  +  «q«/r". 


n) 


(2.141 


i2.1.‘i) 


The  detailed  equations  will  be  presented  in  the  Appendix.  Here  we  emphasize  their  form.  In 
particular,  let  us  regard  q  and  q',  as  parameters.  Then  Eq.  (2.14)  with  Eq.  (2.15)  has  the  form 
(for  the  term  shown  in  Eq.  (2.15)) 


^/I(q^,  t) 

ar 


=  0(qj,  t)  +  K(qi)M(i) 


(2.16) 


where 

\f(t)  =  j H(qj")n{q/'.  Ddq," 


(2.17) 


after  all  parameters  such  as  q,q/  are  suppressed.  Multiplication  of  Eq.  (2.15)  by  Hiq.)  and 
integration  leads  to  the  reduced  ordinary  differential  equation 


=  /(f)  +  AMU) 
dt 


(2.18) 


where 

lU)  =  j Hiq^i^iq^.  Ddq^  ;  A  =  j Hiq^iKiq^^dq,  (2.19) 

2.5.  The  Shape  of  the  Packet 

Not  only  is  the  remaining  computational  task  greatly  simplified.  Eq.  (2.17)  already 
informs  us  that  (aside  from  a  choice  of  nonnalization)  the  phonon  wave-packet  operators  are 

•h 


If  there  are  S  subbands  of  imponance.  there  is  one  packet  operator  for  each  choice  of 
(«',  n)  or  5(.5 -I- 1  )/2  differently  shaped  packets.  For  the  important  case  in  which  only  the 


175 


'SO  M.  Laz  &  It  .  t  'rii 


II  =  0  stales  panicipate  there  is  one  such  packet  (tor  each  q): 

1  I  ,  I  ,s-xp(- 

(/.  I  =  J  I  .III-'  I  ! 

V.  -  'J</-+r/: 


12.21) 


The  shape  ot  this  packet  in  ordinary  share  may  be  obtained  by  multiplication  by  expiii/.ci  and 
mteeratine  over  </.: 

'  J/^l:  I  I-  (2.22) 


\s  here 

Ki: 


.  e\p|i(/.(;  -  r  il  ,  ,i 

:  I  =  j  (/(/.  =  IK  itic/  I  ;  -  :  I  i 

1/: 


I  2,2.^  I 


has  the  torm  of  a  iTKxlitied  Bessel  function.  This  form  arises  from  the  C'oulomb  nature  of  the 
interaction.  If  we  had  used  a  point  interaction  u/'  +</.-)  '  '  would  have  been  replaced  by 
unity  and  /?(r  -  by  6(r  -  :')  so  that  the  packet  shape  would  simply  be  I  l^oC)  I ".  More 
eenerally.  the  packet  shape  is  a  convolution  of  the  electron  density  with 

I :  -  I  )  for  each  q. 

2.6.  Further  Simplifications 

We  have  not  written  the  explicit  form  of  the  phonon-phonon  interaction  by  means 
of  which  the  relevant  longitudinal  optical  phonons  decay  into  acoustic  phonons  because  we 
have  replaced  that  process  by  a  relaxation  process  of  the  form 


(U 


rr 


^lq(r/?*  Ill 

^op 


12.24) 


where  the  decay  time  for  optical  phonons  t.p  has  been  estimated  in  the  experimental  papers''* 
to  be  7  psec. 

For  times  larger  than  a  picosecond  it  has  been  found  by  Monte  C'arlo  calculations"'*  ’’ 
and  by  our  own  quasi-analytical  procedure"'  that  the  electron  distribution  has  equilibrated 
relative  to  two  macroscopic  parameters:  an  electron  temperature  I,  it)  and  a  dnft  veliKity 
v,  i/).  both  of  which  may  be  lime  depeiulent.  By  inirrHlucing  these  parameters,  the  equation 
for  the  electron  distribution  is  replaced  by  ordinary  differential  equations  for  these  parameters. 

I  inally.  since  we  have  been  dealing  with  variables  such  as  n,,iq..q’.)  it  was  con- 
\enieni  to  use  a  wave  packet  construction  on  both  left  and  right  w.ive  sectors,  lluis  our 
reduced  phonon  variables  are 
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n'n.  ^  * 


where 


IL  </.-)  «qW.-.  ‘i  :-  I'  ‘//l 

_ _ _ _ 

^  n'n.  m 


,  ‘f  -M 


>n'm.  n'n^^i  ^  ^  J  i  J  '  '::m  ^ “  1  *Sm^“  1  ^ “  2 


1 2.25) 


(2.26) 


Thus,  m  general,  we  will  gei  coupled  ordinary  dit'terentiai  equations  tor  /',.(n.  v,(/).  and 
V  ■  (/ ) 

M  «.  m  m '  *  '• 

2.7.  The  reduced  Equations 

l  or  simplicity,  wc  shall  write  here  only  the  special  case  in  which  only  ihe  lowest  sub- 
band  contnbutes.  Thus  we  shall  set 


V'otq,  n  =  .VimiHitq.  n 
The  phonon  equation  can  then  be  written 

tWolt/,  I)  .Vo«/.  rt-t-  1  .Vo(</.  O  'l-  V/lF/  ) 


i2,;7) 


(:.2S) 


(II  r,.{t))  T.,p 

where  the  rate  of  phonon  emission  I  /%  and  the  rate  of  phonon  absorption  1/Tj  are  given  by 
1  TO*  1  Fm  (]K)(q)  f^iq) 


Xpiq.  T,.)  Vi  q  leooooiq.bii) 


12  • 


=exp  1  riwikpT^it)]  — — - 

r,.)  i^iq. 

where  T(H)  (Kitqi  was  detined  in  Eq.  (2.26),  and  the  matrix  element  /'*'  is  given  by 

1 


(2.29) 


(2.30) 


7-,(t))=  j  kdk 


q  „  'm<o  Ii.2  _,<7  ?»"i2|l  ■’ 


(i  mw 


1*^  _,X  + 

2  Fiq 


f  {El,,  r,.(n)|l  -  f  {El,  - /iO)f,.  r,.(r))|  . 


(2.31) 


liq.  (2.30)  shows  that  the  ratio  of  rates  of  absorption  and  emission  is  governed  by  the  instan- 
laneous  electron  temperature.  Dynamic  screening  effects  ,ire  contained  in  ihe  ilielectnc 
response  function  foKmlt/.  O)/  )  used  in  Fiq.  (2.29). 
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Fig.  1 .  1/T  as  a  ftjnction  of  ihe  electron  temperature  T,.  for  =0.  2.  5,  and  10  psec  and  (n  =0)  lowest 
subband  occupation.  Dynamic  screening  is  included  (solid  curves).  The  dashed  curve  shows  for  com¬ 
parison  a  calculation  for  •Cop  =  0  when  static  screening  is  used.'^  Comparison  between  the  T=0  (no 
bottleneck  case)  and  t=  10  psec  shows  a  reduction  of  one  order  of  magnitude  in  the  latter  case. 


2.8.  Energy  Loss  Rate  in  Steady  State  Case 

The  first  explicit  calculation  is  made  for  the  steady  state  case.  In  that  case,  the  elec¬ 
tron  temperature  Tg(t)  =  T,  will  be  independent  of  time  and  assumed  given.  The  problem  is 
to  calculate  the  rate  of  energy  transfer  from  the  electron  gas  into  the  LO  phonons,  and  to 
express  the  result  in  the  form  used  by  the  experimentalists.  The  latter  fit  their  experimental 
data  with  an  expression  of  the  form: 


FbT; 


(2..A2) 


Then  1  T  is  plotted  as  a  function  of  electron  temperature. 
We  calculate  1  lx  by  using  the  expression: 
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PeU)  = 


t) 

Ft 


(133) 


Comparison  with  the  previous  equation  yields  1/t.  A  plot  of  l;x  obtained  in  this  manner  is 
shown  in  Fig.  1. 

We  note  that  in  the  case  in  which  T.,p=().  equilibration  of  the  LO  phonon  modes  takes 
place  instantaneously.  Thus  no  bottleneck  effect  will  occur.  But  the  effective  energy  transfer 
rate  is  reduced  by  an  order  of  magnitude  if  one  takes  i=  U)  psec.  The  dashed  curve  demon¬ 
strates  that  neglect  of  screening  would  have  little  effect  on  the  results. 

If  we  apply  the  steady  state  condition: 


tf/Vnlq.  n 

(it 


(2.34) 


we  can  solve  for  Nniq),  the  non-equilibrium  phonon  occupancy  associated  with  the  transverse 
phonon  wave-vector  q. 

Using  the  Planck  formula, 

^'o(q)  =  l/{expi«COt/^Br(q)|- 1}  (2.35) 

the  results  can  be  expressed  in  the  form  of  a  temperature  T  (q)  for  the  phonons  of  a  given 
transverse  wave-vector  q.  The  results  displayed  in  Fig.  2  show  that  the  predominant  heating 
occurs  for  small  wave-vector  phonons. 


Fig. 2  "Phonon  temperature"  Tlq)  as  a  function  of  q  for  .SO  K  </■,.<  3(X)  K  and  lowest  subband 
occupation.'^. 
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We  have  also  made  calculations  when  electrons  and  holes  are  simultaneously 
present,''  and  when  more  than  one  sub-band  is  ix'cupied."'’  A  comparison  between  expen- 
ment  and  theory  tor  this  case  is  shown  in  Fig. 


Fig.  3.  The  carrier  tetnperatune  is  plotted  against  the  energy  loss  rate  per  carrier  for  the  electron  and 
hole  ca.se.  For  the  electron  case,  the  effect  of  hot  phonons  is  found  to  be  an  order  of  magnitude, 
whereas  for  holes  it  is  much  less  than  an  order  of  magnitude.  The  points  arc  the  experimental  data  of 
Shah  ct  al.‘^ 

In  Fig.  3  the  curve  for  electrons  is  obtained  by  including  only  the  polar  interaction  with  LO 
phonons.  The  calculation  for  holes  includes  polar  coupling  to  LO  phonons  and  coupling  to 
both  LO  and  TO  phonons  via  the  deformation  potential.  Only  heavy  holes  are  considered.  It 
is  shown  that  the  hot  phonon  effect  is  strong  for  electron  relaxation,  but  is  weak  for  hole 
relaxation.  This  follows  because  for  holes  the  phonon  emission  rate  1/T>,((/)  for  given  q, 
which  is  similar  to  Eq.  (2.29)  for  \lx,(q),  is  much  smaller  than  that  for  electrons.  From  Eqs. 
(2.29)  and  (2.31),  we  see  that  MXfiq)  is  proportional  to  which  is  proportional  to 

L>«,  with  rn  the  effective  mass.  The  large  effective  mass  of  holes  leads  to  a  weaker  build-up 
of  hot  phonons  by  the  hole  gas.  On  the  other  hand,  heavy  holes  preferentially  emit  phonons 
with  larger  momentum  q.  because  of  the  large  effective  mass  of  hole.  Therefore,  the  phase 
space  of  q  that  contributes  to  the  total  cixiling  rate  for  holes  is  much  larger  than  that  for  elec¬ 
trons.  This  leads  to  the  larger  cooling  rate  for  holes  than  that  for  electrons. 
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2.9.  Time-dependent  relaxation 

(3ur  starting  equations  are  valid  for  the  time-dependent  case.  We  simply  do  not 
assume  time  denvatives  vanish.  The  electron  energy  can  he  wntten  as  a  sum  over  the 
transverse  k  vector: 

<E,U)>  =  ~  I- 

k 


Tor  simplicity,  we  again  consider  here  only  one  sub-band.  The  rate  of  electron  temperature 
change  is  then  given  by: 

^  1  r)<E(t)>  (St?) 

fit  CutT^i/))  i)i 


where  is  the  specific  heat:  C'„  =  d<£(f)>/dT,.(f ).  Since  the  phonon  equations  depend  on 
the  instantaneous  temperature  we  now  have  coupled  ordinary  differential  equations  con¬ 
necting  the  phonon  tK'cupancies  .Volqi  and  the  electron  temperature.  It  is  assumed,  of  course, 
that  r,,(f  =  01  is  given.  A  comparison  is  given  in  Fig.  4  of  our  theoretical  results  with  expen- 
ment 

As  the  electron  gas  relaxes,  its  temperature  decreases  as  shown  by  the  solid  curve  in 
Fig.  5.  The  associated  rise  in  the  phonon  temperature  is  shown  by  the  dashed  curve  which 
merges  with  the  solid  curve  as  the  combined  system  relaxes  to  the  lattice  temperature.  When 
T„p  is  set  equal  to  zero,  the  phonon  bottleneck  effect  disappears,  and  the  electron  temperature 
falls  more  quickly  as  shown  by  the  dot-dash  curve. 


F  ig.  4.  Time  dependent  relaxation  of  an  electron  gas  ol  density  n,  =  Id'*  cm'  starting  at  an  initial  lem- 
(icraiurc  ol  T,{())  =  ^5^)K.  I'lic  optical  phonon  decay  rate  has  been  given  the  accepted  value  ol  t„p  =  7 
pscc.  The  expcnmental  data  arc  from  Ryan  cf  u/  ‘  corresponding  lo  a  .ID  elcctmn  density 
n  =  5'  !()'’'( m  '  and  lo  a  maximum  power  absorbed  by  the  sample  ol  -  .s()  mW. 
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Fig,  5.  A  plot  of  cleciron  icmpcraiure  T,  against  time  (solid  curve  for  T„p  =  7psec  and  dot-dash  curve 
for  =0)  and  of  the  phonon  temperature  Tp(if)  at  <?  =  1.3  x  10* cm'' . 


3.  FIJRTHER  APPLICATIONS  OF  THE  THEORY 
3.1.  Hot  electron  transport 

When  a  strong  dc  field  is  applied,  the  electron  gas  acquires  an  elevated  temperature  as 
well  as  a  drift  velocity.  Moreover,  the  drift  mobility  is  reduced  more  in  the  presence  of  a 
phonon  bottleneck  than  in  its  absence.  This  issue  was  studied  by  Lei  and  Horing^*  and  by 
us‘^  The  methods  described  earlier  involving  the  introduction  of  a  phonon  packet  remain 
valid.  It  is  convenient,  in  addition  to  separate  the  center  of  mass  motion  of  the  electron  gas 
from  its  relative  motion.  The  density  of  electrons  will  again  be  assumed  sufficiently  high  that 
the  electron  gas  can  be  assumed  in  quasi-equilibrium  with  a  given  drift  velocity  and  a  given 
temperature.  A  term,  which  describes  the  center-of-mass  motion  under  an  applied  field  should 
be  added  upon  to  the  Hamiltonian  described  in  Sec.  2.3; 

H,  =  p'  2M  -  NeE  R  .  (3.11 

with  P  and  R  the  center-of-mass  momentum  and  coordinate  operators,  respectively,  N  the 
total  number  of  electrons,  and  M  =  Nm.  with  m  the  electron  effective  mass.  The  electronic 
states  are  labeled  by  the  2D  wave  vector  in  the  relative  coordinate  system,  k,  and  the  discrete 
subband  index,  n;  a„Y  and  (Eiq.  (2.1))  now  are  the  corresponding  electron  creation  and 
annihilation  operators  in  the  relative  coordinates  and  E'„^  is  the  energy  of  an  electron  in  the 
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{n.  k)  state.  Theret'ore,  the  energy  exchange  of  electron  in  the  relative  cixjrdinates  obeys  the 
conservation  condition: 

l^n'k-  -  /-'nk  =  ‘ 

with  the  drift  veUx;ity  of  electrons.  The  intracollisional  effect'"  ",  namely  the  effect  of  the 
electric  field  during  the  course  of  a  collision  can  be  neglected  for  the  fields  considered  in  our 
work'".  If  only  the  lowest  subhand  is  considered,  we  obtain  a  pair  of  time  evolution  equa¬ 
tions  for  the  drift  velocity  v,.  and  the  total  energy  of  electrons  in  the  relative  cixirdinates: 

— r - =  -/V<>E  -  t)  ,  i.v.^) 


d£’. 

(h 


=  -  ( fiO)/ 

o-LO  '» 


q  v,.)£r_LO<q.  >) 


Here 


^e-Lolq.  0  - 


n 


(.V4) 


where  A/o(q,  t)  is  defined  by  Eq.  (2.27)  and 

/Ji/fq,  0)t)  =  Y  XZ/<ik^O  I  1  -/o'kOOl  5^'  SlEok  -  t  liq  v,  +  /icot)  .  1.^.6) 

'''  kk' 

where  the  electron  distribution  function  in  the  relative  coordinates  /ok(0  is  assumed  to  be  a 
Fermi-Dirac  distribution  with  temperature  r«(r).  With  the  above  approximations,  the  problem 
is  again  reduced  to  coupled  ordinary  nonlinear  differential  equations.  Solving  Eqs.  (.T3)  and 
(3.4)  and  corresponding  equation  for  N^,  which  is  similar  to  Eq.  (2.28),  we  can  determine  v,,, 
Tf  as  well  as  A^olq) 

Results  for  the  mobility  reduction  and  temperature  increase  in  a  steady  applied  field 
are  shown  in  Fig.  6. 
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Fig.  6.  The  normalized  mobility  of  electrons.  |i(E)/Ai(E  — >  0)  (solid  curves),  and  the  electron  tempera* 
lure,  Tf,  (dashed  curves)  as  functions  of  the  external  eleetric  held  I E  I  at  =  0  and  T„p  =  7  psee. 
The  crosses  represent  the  experimental  data  from  Fig.  2(c)  of  Keever  et  at  Tl  =  77  K. 


l  ig.  7,  Optical  phonon  icmpcraturc  ",  /"..piq).  .is  a  lunciion  ol  q  lumt:  lo"  cm  '  i  ai 
I  B  I  =  5(X)  V/cm.  The  </,-</.  plane  represenis  T,^  =  77  K.  ihc  peak  value  of  T,^  is  Wh  K.  Insci: 
The  average  wave  vector  of  nonequilibrium  1.(7  phonons.  <  </,>.  as  a  function  of  the  electric  held. 
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Phonon  heatinii  produced  by  the  strong  electric  field  is  not  isotropic,  but  is  more  effec¬ 
tive  for  phonons  whose  propagation  direction  is  in  the  direction  of  motion  of  the  electrons. 
The  phonons  also  acquire  a  mean  ■momentum”  <r/,>.  -A  contour  plot  of  the  phonon  tempera¬ 
ture  rise,  and  an  inset  of  the  mean  momentum  are  shown  in  I'ig.  7. 

We  also  studied  the  tune  evolution  of  electrons  beginning  from  Twitch  on  of  a 
applied  electnc  field  until  arrival  at  a  steady  state.  Fig.  S  displays  the  results.  The  presence  of 
hot  phonons  leads  to  time  delay  in  arriving  at  a  steady  state. 


Fig.  8.  The  drift  velocity  of  elecuons,  Iv,  I  (solid  curvcsl  and  the  electron  temperature.  T,  (dashed 
curves)  as  functioas  of  time  t.  The  parameters  are  the  same  as  those  of  Fig.  h. 

3.2.  Relaxation  of  the  electron-hole  plasma 

When  an  undoped  quantum  well  device  is  illuminated  by  a  high-power  laser,  photo- 
excited  electrons  and  holes  are  produced  simultaneously.  Therefore,  it  is  necessary  to  study 
the  cooling  process  of  a  photo-generated  electron-hole  plasma.  This  process  involves,  not  only 
the  carrier- phonon  coupling,  but  also  the  electron-hole  interaction.  This  issue  was  studied  by 
Pdtz'’  in  bulk  GaAs.  and  then  by  Marchetti  and  Pcitz''*  in  a  GaAs-GaAlAs  quantum  well 
using  the  approach  described  above,  since  hot-phonon  rcabsorption  can  also  slow  down  con¬ 
siderably  the  carrier  cooling  in  the  electron-hole  plasma. 

The  energy  of  photo-excited  electrons  is  assumed  below  that  of  the  I.  valley,  so  inter- 
valley  scattenng  is  not  included.  The  carrier  distributions  are  mtxleled  as  a  time-dependent 
Fermi-Dirac  distribution  functions.  The  temperatures  in  the  distnbution  functions  for  electrons 
and  holes  are  allowed  to  be  different  to  account  for  noninstantaneous  energy  transfer  between 
two  carrier  systems.  The  equations  for  the  time  evolution  of  total  electron  energy  tE,.i  and 
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hole  energy  {£/,)  is  derived  as 


f)£(. 

IT 


(3.7) 


and 


IT 


+  Rf.t,-Rt,^lQ-R),^jQ  +  R 


h 

L 


(3.8) 


Here  Rf.),  is  the  power  loss  by  electrons  to  holes  via  the  screened  Coulomb  interaction. 
/?,._L()  and  /?/,_Lo  are  the  rate  of  energy  loss  by  carriers  to  LO  phonons  via  the  polar  Frohlich 
coupling,  Rt,-To  is  the  rate  of  energy  loss  by  hole  to  TO  phonons  via  the  deformation  poten¬ 
tial  coupling.  Rl  and  /?*  are  the  power  input  from  the  laser  to  electrons  and  holes,  respec¬ 
tively.  Only  the  lowest  subband  in  the  quantum  well  is  considered. 

The  corresponding  time  evolution  equations  for  non-equilibrium  LO  and  TO  phonons 
(similar  to  Eq.  (2.28))  are  coupled  to  Eqs.  (3.7)  and  (3.8).  Therefore,  the  temperatures  T,.  T* 
and  the  phonon  distributions.  N(f^(q)  and  can  be  determined  as  function  of  time  t. 
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Fig.  9.  Electron  and  hole  temperatures  as  functions  of  time  for  n;n=0.5xl()'‘ cm‘‘  and  four  values  of 
well  thickness:  £  =  25 A.  solid  line;  £  =  5()A.  dot-dashed  line;  £  =  UX)A.  da.shed  line;  £  =  200A,  dotted 
line.  The  four  curves  starting  at  r  =  2(X)K  represent  hole  temperatures,  the  others  represent  electron 
temperatures.  This  figure  is  Fig.  1  from  Marchetti  and  Pdlz.” 
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In  Fig.  9,  the  electron  and  hole  temperatures  are  displayed  as  functions  of  time  for  dif¬ 
ferent  layer  thickness  L.  The  exchange  of  energy  between  electrons  and  holes  plays  an  impor¬ 
tant  role  in  the  initial  stage  of  the  relaxation.  The  main  ponion  of  the  photon  excess  energy  is 
given  to  the  phoiogenerated  electrons.  Initially,  the  kinetic  energy  of  excited  holes  is  below 
the  threshold  for  optical  phonon  emission.  However,  the  e-h  coupling  rapidly  transfers  energy 
from  electrons  to  holes  and  thus  ensures  significant  participation  of  the  holes  in  the  ctxiling 
process.  This  energy  loss  channel  is  also  imponant  after  the  buildup  of  LO-phonon  modes 
which  couple  to  electrons  slows  the  cooling  of  electrons.  Reabsorption  of  phonons  becomes 
important  within  less  than  a  picosecond  after  the  onset  of  the  laser  pulse  and  leads  to  a  con¬ 
siderable  reduction  of  the  earner  cooling  rate  in  the  later  stage  of  the  relaxation. 

One  can  see,  from  Fig.  9,  that  the  carrier  cooling  rate  is  rather  insensitive  to  vanations 
in  the  well  thickness,  if  the  sheet  density  is  kept  constant.  This  result  is  in  agreement  with 
experiment. .  On  the  other  hand,  their  calculation^**  indicates  that  for  given  layer  thickness 
and  times  beyond  0.5  ps,  cooling  occurs  at  a  slow  rate  at  higher  values  of  sheet  density, 
because  of  a  strong  buildup  of  nonequilibrium  optical  phonons  at  higher  sheet  density. 

3J.  Electron-hole  Transport  and  Negative  Mobility 
The  luminescence  measurements  of  the  photoexcited  electron-hole  (e-h)  plasma  in 
quasi-2D  quantum  wells  and  the  relative  theoretical  study,  which  was  discussed  in  Sec.  .T2. 
are  important  to  derive  information  about  relaxation  of  this  system.  On  the  other  hand,  tran¬ 
sport  measurements  of  a  photoexcited  electron-hole  system  under  an  applied  electric  held 
have  also  provided  some  interesting  results.^*’^’  Recently,  Hopfel,  Shah.  Wolff,  and  Gos- 
sard^^  found  that  in  a  such  system  the  minority  electrons,  which  are  injected  by  laser  pumps 
on  the  p-modulation-doped  quantum  wells,  can  move  in  the  direction  of  the  external  electnc 
held.  This  negative  absolute  mobility  of  electrons  occurs  because  of  strong  electron-hole 
drag.  This  subject  was  theoretically  studied  by  Cui,  Lei,  and  Horing^*  and  by  us*^. 

We  first  briefly  discuss  the  condition  for  negative  absolute  mobility  of  electrons  in  the 
region  of  weak  electric  field  where  the  conductivity  is  linear.  In  the  steady  state  we  have  the 
force  balance  equation  for  carriers: 

-  P‘-'' -  =  0  .  (  .T9  ) 

where  E  is  the  external  electric  field.  represents  the  frictional  force  due  to  the  earner- 
lattice  interaction  and  P‘“''  represents  the  frictional  force  upon  the  p-type  earners  due  to  the 
carrier-carrier  interaction  with  the  v-type  carriers.  It  is  obvious  that  F"**  ''  =  -  F"'  In  the 
region  of  linear  conductivity  we  have 
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(4. 10) 

=  Uu«v4^'''(Vu  -  Vy)  , 

(4.11) 

where  is  the  drift  velocity  of  jj.  type  of  carriers,  ,4^'*'  represents  the  contribution  to  the 
resistivity  (per  earner)  from  p-l.  scattenng,  and  A**  ''  relates  to  the  contribution  to  the  resis¬ 
tivity  (per  carrier  (i)  from  p-v  scattering  normalized  to  per  carrier  v.  From  Fqs.  (44))  -  i  4.1 1) 
we  immediately  obtain  the  mobility  for  electrons,  p^,: 

(n,^.  —  /i^)4*'  ^  +  4^  ^ 

fio  =  — i — iri - — i - ir“in-  (4.12) 


[iquation  (4.12)  indicates  that  for  electrons  a  negative  absolute  mobility  is  only  possible  when 
<  «h.  (The  corresponding  statment  for  hole  mobility  could  require  /i|,  <  /i^.i  It  is  deter¬ 
mined  by  the  difference  between  the  density  of  electrons  and  the  density  of  holes,  and  the 
competition  between  electron-hole  drag  and  hole-lattice  scattering.  At  low  temperature  and 
under  a  weak  electric  Held,  the  former  dominates,  so  the  mobility  of  electrons  is  negative. 
When  the  lattice  temperature  or  the  electric  Held  increase,  the  latter  tends  to  dominate  and  the 
mobility  of  minority  electrons  becomes  positive. 

The  carrier  dynamics  is  derived  in  a  similar  way  as  discussed  in  Sec.  4.1  and  Sec.  3.2. 
We  separate  the  center-of-mass  moHon  of  each  type  of  carrier  from  its  relative  motion.  The 
2D  momentum  for  p-type  carriers  in  the  relative  coordinates  is  deHned  as  =  Sk  - 
and  a  Fermi-Dirac  distribution  function  is  assumed  at  temperature  T^^.  in  the  relauve 

coordinates.  The  exchange  of  momentum  and  energy  in  the  relarive  coordinates  can  then  be 
written  as 

k  -  k'  =  q  ,  -  E^\-  =  K(  <a  -  q-v^,  )  .  (4.14) 


A  set  of  coupled  equations  for  the  time  evolution  of  the  drift  velocities  of  the  center-of-mass 
and  energies  in  the  relative  coordinates  for  each  type  of  carrier  can  be  derived.  We  obtain 


dr 


L 


(4.14) 


and 

-^i-L  -^i-v 

<)E  (f)  _o£  (n  <)E  (/) 

- 5r~ 

in  Eq.  (4.14)  represents  the  frictional  force  due  to  carrier- lattice  (impurity)  interaction. 
The  expression  of,  for  example.  LO  phonons  is  similar  to  Eqs.  (4.4)  -  (4.6).  Mere,  the  earner 
interactions  with  LO,  TO.  acoustic  phonons,  and  impuniies  are  considered.  .Also,  the  hot- 
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phonon  effects  for  1.0  and  TO  phonons  are  included.  F**  which  represents  the  fnctional 
force  due  to  e-h  scattering,  with  v  a  different  type  of  carrier  than  p.  is  given  by 


Ik 

h 


III  I  I  '  i  ^  o'*!-  I  qVp-t 

q  i'.i  I  k'.k  p'.p 


-dl/ffk  -t-  Tip  -  /Vk-  -  l  ,p  +  /iq  (v,,-v.,)i  . 


I lb) 


where  the  term  /iq  (v^,-Vv)  appears  in  the  last  6  function  because  the  energies  T  for  electrons 
and  holes  are  defined  in  different  coordinates,  and  (o  =  -  /-Vk  •  ^•*1' 

1  V  1.^1  represents  the  energy  loss  rate  in  the  relative  system  due  to  e  h  scattering.  .An  equation 
for  the  rate  of  energy  loss  can  be  obtained  by  replacing  in  Hq.  (.Tib)  by  -oT  if)  di 
and  the  prefactor  /iq  on  the  right  hand  side  of  Kq.  i.’.lb)  by 

-  4V 

The  term  I  ,  ,(q.  w)  in  Kq.  i.vlb)  represents  the  2D  screened  scattering  iiiatnx 
between  earners.  Here,  one  should  be  careful  in  treatment  of  screening  effect  when  two  kinds 
of  earners  coexist.  The  unscreened  electron-electron,  hole-hole,  and  electron-hole  scattenngs 
occur  via  Coulomb  interactions  given  by 


27Wpev 

EnqA 


(.T17) 


where  the  form  factor  y'y**/*  given  by  Eq.  (2.2b).  and  en  is  the  stattc  dielectric  constant. 
.Since  the  strength  of  the  electron-hole  interaction  is  of  the  same  order  as  the  elecn-on-electron 
interaction  and  the  hole-hole  interactions,  the  dynamic  screening  effect  should  be  expressed 
by  a  more  complex  form  than  that  in  the  case  of  a  single  type  of  carrier.  These  interactions 
are  shown  in  Fig.  10  in  the  random  phase  approximation  (RP.A).  ,\ccording  to  F'ig.  10.  the 

~  HV 

screened  carrier-carrier  scattenng  matrix.  V'  .  satisfies  the  following  equation: 

(■bv  ^  ^.pv  ^  ,  i3.is) 

n 


where  H'^  is  a  diagonal  matrix,  its  elements  are  the  density-density  correlation  functions  for 
2D  earners.  If  we  define 


■■  oc 

*  oh 

r 

1 

1 

11  = 

1 F  0 

[v""  l'"'' 

! 
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rile  screened  camer  camer  scattering  matrix  can  then  be  expressed  in  matrix  tomi: 


T  =  I  1  -  I  11  I  '  I 


i.v20i 
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When  only  the  lowest  hand  tor  electrons  and  holes  are  occupied.  I  can  be  straightlorwardly 
written  as 


t  ig.  10  Diagrams  lor  scmening  etiect  ol  eamer-camer  potential  in  the  RPA.  The  single  and  double 
hon/ontal  solid  lines  represent,  respectively,  the  unscreened  and  screened  camer-carrier  intcracuon. 
V'**''  and  V'*”"  The  bubble  represents  the  density-density  correlation  function.  IT'. 

Solving  a  set  of  coupled  equations.  Eqs.  (3.14)  and  (3.15)  and  corresponding  equations 
for  hot-phonons,  we  can  determine  the  v,,  v*.  T,,  as  well  as  the  distnbution  of  hot-phonon 
wavepackets  (LO.  TO).  The  mobilities  for  both  earners  in  a  weak  electric  field  are  shown  in 
l  ig.  1 1  as  functions  of  lattice  temperature,  T^.  Fig.  12  shows  the  weak  field  mobility  of  elec¬ 
trons  and  holes  as  functions  of  the  electron  density. 
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Fig.  1 1 .  Mobility  of  electrons,  itnd  mobility  of  holes  .  gi,-  •i''  lunciions  ol  lattice  lempcralure.  F, ,  m 
a  weak  electric  held.  Data  tor  electron  mobility  come  from  Hoplel  et  al'"  . 


Fig  12  Mobility  ol  electrons,  g^  (solid  curvcsi.  and  mobility  ol  boles,  ulashcd  curvesi,  as  lunc 
lions  ol  dcnsiiy  ol  electrons.  ai  dilletvnt  lattice  temperaturc.  If  i  in  a  weak  electric  Held  i 
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4.  OTHER  theoretical  WORKS 


4.1.  Monte  (.'uriu  Simulation  of  N«tnequilibrium  Phonon  etTect 

In  ihe  above  discussion  a  simplifying  assumption  for  the  carrier  liistribuiion  function, 
namely,  that  is  a  l  errni-Dirac  function  at  the  earner  temperature  T.  is  assumed.  This  assump¬ 
tion  usually  is  a  good  one  for  the  high  earner  density  and  for  the  time  scale  larger  than 
pecosecond.  However,  recent  progress  in  ultrafasi  expenments  leads  to  processes  in  which 
the  earner  distnbution  can  not  be  regard  as  in  i|uasi  equilibrium.  I'se  of  Monte  ('arlo  tech¬ 
niques  to  solve  the  dynamical  equation  on  the  supercomputer  is  the  way  of  obtaining  the 
nonequilibrium  earner  distribution.  Lugli  ft  ill  have  presented  a  senes  ibeoreiical  results  on 
nonequilibrium  phonon  effects  based  on  a  novel  Monte  Carlo  algorithm.'"  They  have 

calculated  the  results  for  both  the  hulk  and  quantum  well  cases.  In  the  case  of  .i  quantum 
well,  the  effect  of  several  subbands  has  been  included.  However,  a  simple  intxiel  vvas  used. 


that  is,  that  the  t/j,  components  of  the  phonon  distribution  are  localized  in  a  region  of  wave- 
vector  space  of  extent  1  /L.  with  L  the  width  of  well. 

The  physics  of  the  dynamical  evolution  of  the  carrier-phonon  system  is  very  similar  to 
what  we  have  discu.ssed.  The  coupled  Boltzmann  equations  are  given  by 


^>/k 

dt 


|4.1) 


rw;/’ 

fJt 


(4.2) 


where  the  superscripts  i  and  /  indicate,  respectively,  the  type  of  earner  (electron  or  holes)  and 
of  phonon  modes  (LO,  TO,  .  .  .)  considered.  {dN^)ldt  I  ph-ph  's  obtained  by  relaxation  time 
approximation,  given  by  Eq.  (2.24).  Fig.  1.^  shows  the  evolution  of  the  electron  total  energy 
as  a  function  of  time  in  an  n-type  quantum  well  during  and  after  the  laser  pulse.  The  excited 
electrons  lose  energy  mainly  through  the  interaction  with  the  background  electrons  and 
through  the  emission  of  FO  phonons.  A  much  slower  relaxation  is  found  when  nonequili¬ 
brium  phonon  effects  are  included.  Fig.  14  shows  the  LO-phonon  distnbution  .it  different 
times. 
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(•ig.  LV  Average  electron  energy,  measured  in  equivalent  temperature,  as  a  lunction  of  lime  dunng  and 
after  the  la.ser  excitation  lor  two  different  excitation  sheet  densities.  The  position  of  the  energy  levels  m 
the  well  is  shown  in  the  insert.  This  ligure  is  Fig.  12  from  Lugli  el  al.^‘ 


l  ig.  14.  LO-phonon  distnbuiion  as  a  lunction  ol  total  parallel  momentum  ilor  i/.  =ih  lor  times  during 
and  after  the  laser  excitation.  This  ligure  is  Fig.  14  from  Lugli  et  al.'' 
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When  ihe  laser  photon  energy  is  about  2  eV.  photon-exeitation  in  bulk  CiaAs  includes 
iransition  trom  three  hole  bands.  The  energy  distribution  of  phoio-e.xcited  electrons  is  charac- 
tcri/ed  by  three  distinct  peaks,  as  shown  in  the  insert  of  l  ig.  l.'s.  In  liaAs.  about  bO'f  of  the 
photiK’.xcited  carriers  transfer  to  the  satellite  valleys  dunng  the  laser  pulse  (the  average  time 
for  1'-/  tt.insition  via  phonon  emission  or  absorption  is  about  SO  fs).  Camers  return  slowly 
lo  the  r  valley,  with  charactenstic  time  of  2  or  ps.  because  of  smaller  effective  mass  in  the 
r  valley.  I'his  leads  a  slow  rise  of  luminescence  in  CiaAs,  .is  shown  in  1-ig.  1,S.‘'  In  l  ig.  I.S.  a 
comparison  to  the  InP  case  is  also  provided,  where  the  1.  valley  are  kx.aied  at  a  much  high 
energy,  ,ind  do  not  signiticantly  contribute  to  the  ctH>ling  process.  Bv  litting  data  with  the 
Monte  Carlo  calculation,  it  is  detemiined  that  the  I'-/,  defonnation  potential  Oyi,.  is 
1(1.,^  I  lO'^  cV/cm.  'I'his  P-/,  exchange  depletes  the  P  electron  distribution  in  Cia/\s  in  the 
high-energy  region,  above  0..^  eV,  as  shown  in  Pig.  I(i. 
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l  ig.  1.5  l.uminescence  intensity  at  '(X)  K  vs  delay  la)  CiaAs  at  1  4.5  eVi  lexpenmentsi  (hi  the  solid 
curve  IS  ihc  same  as  lai  on  expanded  scale;  the  tilled  circles  show  the  spectrally  integrated  intensity; 
ici  InP,  same  as  (b);  Results  ol  ensemble  Monie  Carlo  calculations  are  also  shown:  ib) 
/)i  1  =  i4,  h,  and 8 1 V  lO"  eV/cm  lor  CiaAs  ilo/enges,  eros,ses,  and  open  circles,  respectively!  and  (c)  InP. 
I'his  ligure  is  Pig.  2  from  Shah  ct  ul  " 
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l  ig,  16,  riic  Monic  Carlo  results  ot  the  Energy  distnbution  tunciion  lor  I -valley  eleeirons  m  la)  GaAs 
and  (b)  InP  at  three  different  time  delays  after  the  excitation.  The  excitation  density  is  5vl()‘''cm  ‘  and 
the  lattice  temperature  is  ,1(X)  K.  This  ligurc  is  Fig.  16  from  Lugli  et  al.*' 


4.2.  Effect  of  LO-phonon  renormalization 

A  series  of  works  by  Das  Sarma  and  his  coworkers'*^"^*’  studied  the  hot-electron  relax¬ 
ation  by  emission  of  LO  phonons  in  both  2D  and  ,^D  cases.  In  addition  to  the  hot- phonon 
effect  and  dynamical  screening,  they  found  that  many-body  renormalization  of  the  LO  pho¬ 
nons  also  plays  a  crucial  role  in  the  power-loss  process  at  low  elecffon  temperatures.  Due  to 
phonon-plasmon  coupling,  the  density  states  of  an  LO  phonon  has  three  branches:  ( 1 )  the 
bare-phonon-like  branch  ne-ar  the  baie-phonon  energy.  t2)  the  pla.smonlike  branch  near  the 
plasmon  energy,  and  (3)  the  low-energy  quasiparticle-excitation-like  branch  in  the 
quasiparticle-excitation  region.  Even  though  the  oscillator  strength  of  plasmonlike  and  the 
quasiparticle-excitation-like  branchs  are  extremely  small,  they  dominate  the  power-loss  pro¬ 
cess  at  low  enough  electron  temperatures  that  most  electrons  have  energy  below  the  threshold 
for  bare  LO  phonon  emission.  This  produces  an  enhancement  of  the  power  loss  at  low  elec¬ 
tron  temperatures  by  many  orders  of  magnitude  relative  to  the  power  loss  to  hare  LO  pho¬ 
nons. 

The  phonon  density  of  states.  A  {q.  (O)  is  given  by 

(0)  =  -it  'imOti/. W)  i4,.'i 

where  D  is  the  phonon  propagator.  For  bare  phonons.  D  is  replaced  by  D'\ 
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/V’l(/.  coi  = 


-M/. 
0)-  -  0)/ 


1 4,4 1 


.ind  liii/4"((/.ci))  =  -7t|6(co-(i)/  )-6(0)  +  (i)/ 41.  which  has  vanishinc  wc-ight  everywhere  except 
■It  the  hare  energy.  By  introducing  the  phoiKin-plasmon  coupling  in  the  random  phase  approxi¬ 
mation  (RP.'\).  the  renonnalized  phonon  propagator  is  obtained  by 

2  CO/ 

Dll/.  COI  =  , - - ^ — , - .  i4..x) 

CO"  -  CO/  -  2o)/  iV/;,  I  !(</.  co> 


where  .V/,,  is  the  electron-phonon  m.itrix  clement,  and  lilt/,  coi  is  the  electron  density-densiiv 
correlation  tunction.  F  ig.  17  displays  the  renonnali/ed  phonon  spectral  tunction  tor  the  .’D 
and  2D  cases.  Three  branches  ot  phonon  nnxies  are  shown  ihai  the  scale  Tacior  y  in 

i  ig,  17  is  different  for  different  curves.  The  strength  >)f  the  plasmonlike  phonon  is  very 
small,  the  t|uasiparticle-e.xcitation-like  iQI'L-like)  mode  is  also  weak.  However.  Ii  is  found 
that  the  QPE-like  phonons  dominate  the  hot-electron  power-loss  prex'ess  at  low  electron  tem¬ 
peratures. 


ut(meV) 


Pig,  17.  (a)  Phonon  spectral  function  for  three  wave  vectors  in  a  .VD  electron  gas  at  an  electron  tem¬ 
perature  of  T  =  .‘iOK  (T/  =0).  For  (>  =2.2x10' cm  '.  the  QPE  branch  exicnds  from  0  to  7  meV.  and 
ihere  are  6-lunction  peaks  near  the  plasmon  energy  14  meVi  and  ihe  phonon  energy  i  'b  X  mcVi, 
.vtiidi  .ire  run  shown  explicitly  in  order  lo  avoid  loo  much  vlelail.  l  or  c4=vri-|(Tem  die  i2PF 
branch  extends  lo  14  meV.  and  there  is  again  a  6-lunction  peak  inoi  showni  at  (w  For 
rj  =  I  ’  ■  10'’ cm  '  the  plasmonlike  phonon  is  l.andau  dam(X'd.  .ind  Ihe  DPF-hke  phonons  exiend  all 
wav  up  lo  (I),  The  electron  density  is  It)'  cm  ’.  giving  a  Pemii  vector  ol  I  4-10'’cm  '  ibi  Phonon 
spectral  lunction  lor  iwo  diHcreni  wave  vector  in  a  2D  electron  gas.  l  or  i/=ii2  7-  lo'em  ‘  ihere  is  a 
(i  liinclioi)  peak  inol  shown)  at  lO/  The  electron  density  is  5- 10'' cm  "  giving  e  Pemii  vector  ol 
I  s-  Mf'cm  ‘  This  Figure  is  Fig  4  Irom  Das  Sarma."’ 
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At  low  temperature  (about  50  K)  the  number  of  electrons  which  can  emit  phonon  with 
bare-mtxie  phonons  becomes  extremely  small.  However,  phonon-plasmon  coupling  creates 
new  low-energy  modes,  even  though  its  strength  is  rather  small,  they  could  make  important 
contribution  in  energy  loss  at  low  electron  temperature.  In  this  case.  Hq.  (2..f2)  tor  the  energy 
loss  rate  is  no  longer  valid  since  it  was  derived  by  considenng  only  bare  1.0  phonon  iiuxles. 
Instead,  Das  Sanna  at  al  derived  the  following  expression  tor  power  loss  by  generalizing  the 
Kogan  formula*’  to  including  phonon  self-energy  correction: 

=  — MiV/‘|n/^,<,„-n/(a))|lmn(</,  0))lmD(t/,  (0)  ,  i-t.6) 

where  «/  =  |exp(co//:fl7t-  1 1  ' . 

In  order  to  include  the  hot-phonon  effect,  one  must  first  rewrite  hq.  i4.(ii  as 

=  5^  J  d(i)a)/?(t/,  co)«7«i))  , 

then  replace  Riq.ta)  in  the  integrand  by  /?(</,  w)/|l -t-T„p/?(t/, coil.  Das  Sarma  <’i  al  suggest 
that  only  the  bare  phonon  mode  is  important  for  hot-phonon  effects. 


I  tg  IS.  Power  loss  as  a  function  ol  the  inverse  electron  temperature  lor  sD  electmn  gas  (solid  linesi. 
ihe  dashed  lines  correspond  to  power  loss  to  bare  l.O  phonons  only.  Ihe  doted  lines  give  the  power 
loss  lo  acoustic-phonon  modes,  l  or  dotted  curves,  the  upixmiosl  curve  iv  lor  111'' cm  .  Ihe  middle 
one  lor  111''’ cm  and  the  lowest  one  lor  Kl'  -cm  *.  T;  =».  and  the  power  loss  is  expressed  in  watts 
per  earner  Iliis  ligurc  is  fig  .s  Irom  Das  Sarnia  <  i 
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Fig.  19.  Power  loss  as  function  of  inverse  electron  temperature  for  two  quantum  wells  (width  SO  and 
2(X)A)  calculated  in  3D  phonon  approximation  (3DPA)  aixl  in  slab-model.  ri,=0,  and  the  power  loss  is 
expressed  in  watts  per  carrier.  This  figure  is  Fig.  1  from  Das  Saima  et  al.*^ 

At  high  temperatures  (T>50  K)  the  bare  LO  phonons  dominate  the  power  loss,  at 
intermediate  temperatures  {20 <T S50  K)  power  is  lost  predominantly  to  the  plasmonlike  or 
QPE-like  branch  of  the  LO-phonon  spectrum,  and  at  still  lower  temperatures  {T  <20  K)  the 
acoustic  phonons  become  important.  Fig.  18  shows  the  power  loss  for  3D  electron  gas  at  lat¬ 
tice  temperature  Ti=0.  Fig.  19  shows  the  results  in  quantum  wells. 

Recently,  M.  W.  C.  Dharma-wardana'**  deals  with  the  same  subject  using  nonequili¬ 
brium  Green’s  function  approach.  He  considers  not  only  the  coupled-mode  spectral  density, 
but  also  the  coupled-distribution  functions.  In  contrast  to  Das  Sarma  and  co-workers,  he  found 
that  the  energy-loss  rate  is  significantly  suppressed  by  coupled-mode  formation  and  quasiparti¬ 
clelike  modes  do  not  contribute  to  the  energy-loss  rate.  The  underlying  physical  picture  is  that 
although  the  coupled-mode  phonon  spectral  density  has  some  weight  in  the  quasiparticlelike 
region  of  energies,  these  coupled-phonon  modes  are  effectively  at  electron  temperature,  not  at 
lattice  temperature.  Hence  there  is  no  enhancement.  Nevertheless,  coupled-mode  phonons  are 
hotter  than  bare  phonons  and  this  leads  to  a  suppression  of  energy-lose  rate. 
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5.  SUMMARY 

We  have  reviewed  in  this  paper  the  el't'eci  of  reabsorption  of  noncquilibnuin  optical 
phonons  on  the  carrier  relaxation  and  transport.  Both  experiments  and  theoretical  study 
continm  that  the  build-up  of  hot-phonons  is  mainly  responsible  for  reducing  the  cixiling  rate 
of  photo-excited  electrons,  while  the  dimensionality  and  dynamical  screening  effect  play  only 
a  minor  role.  In  the  case  of  heterojunctions  or  quantum  wells,  our  approach  can  detenmne 
the  accurate  shape  of  the  phonon  wave  packet  and  produce  a  simple  priKcdure  for  calculating 
the  effect  of  nonequilibrium  phonons.  Our  results  are  in  quantitative  agreement  with  expen- 
ments.  On  the  other  hand,  the  ciwling  of  photo-excited  earners  is  a  nch  field,  which  may  in¬ 
volve  the  electron-hole  (three  kind  of  holes)  interaction,  multi-valleys  and  multi-subband  tran¬ 
sitions,  multi-scattenng  processes,  from  very  high  temperature  down  to  very  low-temperature. 
As  examples,  we  refer  to  the  effects  of  V-L  transitions,  which  delay  the  build  up  of  hot  I' 
electrons  in  a  ultrafast  process,  and  what  is  the  role  of  phonon-plasmon  modes.  There  are 
many  other  factors,  which  may  affect  the  relaxation  process,  such  as  phonon  modes  (interface 
and  slab  modes  in  quantum  wells)^'  and  LO  phonon  relaxation  mechanism. This  subject  is 
worthy  of  further  study. 

ACKNOWLEDGMENT 

The  work  in  the  City  College  of  City  University  of  New  York  was  panly  supported  by 
the  U.  S.  Army  Research  Office  and  the  U.  S.  Department  of  Energy. 

Appendix 

In  this  appendix  the  complete  expressions  for  Eqs.  (2.13)  and  (2.14)  are  given  as  fol¬ 
lows: 

<1 

U  /n  k'IOl  1  -/nk(')lSk'.k«iS(£fiV 

“./nkl  t  ll  i  “/n'k'ft  )|5k'.k-<t^^*nk  “  ^n'k'  “  fltO^  )  ^  .t)| 
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